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ABSTRACT 


The results of numerical calculations of the general-instability strength 
of ring-stiffened circular cylinders are presented in graphical form. The calcu- 
lations are based on Kendrick’s ‘‘second solution’’ which is published in Naval 
Construction Research Establishinent Report No. 244 (Part III). The collapse 
pressures from these graphs agree within 10 percent with those computed by 


Kendrick’s theory throughout the normal range of submarine geometries. 


INTRODUCTION 


One consideration in the structural design of submarine pressure hulls is the possi- 
bility of collapse by general instability, i.e., large deformations of frames and shell between 
holding bulkheads. It is also recognized that the elastic general-instability pressure governs 
the extent to which imperfections reduce the load-carrying capacity of the frames. For this 
reason, an accurate determination of this pressure can be useful in frame design even though 
it is far greater than the pressure encountered at normal operating depth. 

The most reliable of several theoretical investigations of the general-instability prob- 
lem appears to be the recent work of Kendrick conducted at the Naval Construction Research 
Establishment, Rosyth, Scotland. This was published in the form of two separate analyses.) 
Both employ the same general approach, but Reference 2 (Kendrick, Part III) is a more general 
treatment than Reference 1 (Kendrick, Part I) and always gives a lower collapse pressure. 

While Kendrick’s analysis is a significant advance in the study of stiffened cylinders, 
its application to submarine design is difficult because of the lengthy calculations required. 
Recently Bryant,? also of the Naval Construction Research Establishment, developed an 
approximation which agrees closely with Kendrick’s Part I solution. Because of its sim- 
plicity, Bryant’s formula is a valuable aid in design calculations. 

These solutions have been examined at the Taylor Model Basin, and several tests* 
have been conducted with machined models to provide experimental evaluation. In general, 
agreement between experiment and theory was good, the failure pressures usually being slight- 
ly lower than the prediction of Kendrick Part I or Bryant, and slightly higher than that of 
Kendrick Part III.* The difference in the pressures given by these solutions is insignificant 
for cylinders with light frames but becomes much larger for heavy frames such as used on sub- 


marines. It was found, in fact, that for contemporary submarine geometries, Bryant’s formula 


1 References are listed on page 15. 


* It should be noted that this analysis contains two solutions. The first represents a physically impossible 
buckling configuration but is presented because it gives a lower collapse pressure. In this report, all mention 
of Kendrick Part III will be confined to the second solution of that analysis. 


and Kendrick’s Part I solution give pressures as much as 35 percent higher than the Part III 
solution. In such cases use of the more exact solution in design work is desirable, but it is 
impractical because of the extensive calculations required. 

In view of these difficulties, it seemed worthwhile to look for a short method of approx- 
imating Kendrick’s Part III solution. Accordingly, an extensive program was begun at the 
Model Basin to obtain numerical solutions over a wide range of geometries, the objective be- 
ing to summarize the results in some graphical form which would be of practical use to the de- 
signer. Such a presentation would not only provide a quicker and more accurate means of 
determining the general-instability strength of a structure but would present a better picture 
of how variations in the scantlings affect this strength. Because of the very large amount of 
computation involved, this program would have been virtually impossible without the aid of 
the high-speed computer UNIVAC.* 

The results of the calculations are summarized in this report in the form of graphs 
which relate general-instability strength to variations in frame size, frame spacing, shell 
radius and thickness, and compartment length. All calculations were for externally framed 
steel cylinders with a Young’s modulus of 30 x 10° psi. Since the elastic general-instability 
pressure is directly proportional to the modulus, these results are readily applicable to other 
elastic materials having Poisson’s ratio y = 0.3. Moreover, ‘since internal frames theoreti- 
cally provide slightly higher general-instability strength than external frames of the same 
dimensions, the results can be safely applied to internally framed cylinders. The accuracy 
of the graphical results is demonstrated by a comparison with numerical solutions of Ken- 
drick’s Part III theory for a wide range of geometry. The use of the graphs is illustrated in 
Appendix A by a numerical example. In Appendix B, approximate formulas are given whereby 
a frame strength parameter can be determined, and several numerical examples are provided 


to demonstrate the accuracy of the formulas. 


METHOD 


Since Kendrick’s Part III analysis. cannot be reduced to a simple algebraic expression, 
the problem of presenting the theory in graphical form must be approached somewhat indirectly. 
The method adopted consisted of plotting the results of many numerical calculations against 
various geometrical parameters until, by a process of trial and error, a system of coordinates 
was found in which the points followed closely a set of single-valued curves. 

In attempting to define general-instability strength in terms of the shape and size of 
the structure, at least five quantities must be considered, i.e., shell radius, shell thickness, 
compartment length, frame spacing, and some measure of frame size. This situation is further 


complicated by the variability of frame shape and the fact that the number of circumferential 


* Numerical results for more than 200 different geometries were obtained on UNIVAC through the solution 
of a fifth-order matrix by iteration methods. 


bo 


waves into which the cylinder buckles (also a function of the relative dimensions) is a vital 
factor in the determination of general-instability strength. Thus it is clear that the problem 
would be greatly simplified if a few fundamental parameters could be obtained by combining 
some of the dimensions in a rational manner. 

It is shown in Bryant’s paper? that the general-instability pressure can be treated 
approximately as the sum of two terms, one involving the strength of the shell and the other 
the strength of frame per unit length of shell. This latter term was defined as the moment of 
inertia about the centroid of a section comprising one frame plus a length of shell equal to 
one frame spacing. Numerical results indicate that Kendrick’s second solution in Reference 2 
can be broken down fairly successfully in this way. However, a slightly different parameter 
used by Bijlaard? in treating the same problem lent itself better to a graphical presentation and 
gave less scatter in the results. Bijlaard expresses frame strength as the moment of inertia /, 
about the centroid of the combined section of a frame plus an effective length L, of shell. L, 
is taken to be 1.57 RA so long as the frame spacing L, exceeds 2 VRh. For smaller values 


of Ly, L., can be determined from Table 46 of Reference 6. The quantity /, can be written 


e 


A,e2 L.h3 
[ee IT Be (1 
e A if 12 
eal 
ibe 


Where Ff is the radius to the median surface of the shell, 
h is the shell thickness, 
Ay is the frame area, 
If is the moment of inertia of the frame about its own centroid, and 
e is the distance from the median surface of the shell to the centroid of the frame. 
With the parameters as defined above, the general-instability pressure p,, of a stiffened 


cylinder can be expressed as the sum of two terms: 


a 2 h n (2 n [2] 
Pos = 18, Rp’ PR’ Jey 1” ) 


where n_ is the number of circumferential waves, 
L, is the bulkhead spacing, and 


p, is a linear function of h/R. 

re 
100A 
L,/F for values of n of 2, 3,4, and 5. Similarly, Figure 2 shows the variation of Pf with 


The quantity p, can be determined readily from Figure 1 in which p,- is plotted against 


lf, /L,R* for the different values of n. These curves were drawn to fit a large number of 


calculated points.* The same information is presented more concisely in Figure 3 where the 


results from Figures 1 and 2 are combined in one graph having the coordinates 


*Since only values for Pop were obtained from the calculations, p,, was plotted against 1,/L Ro, and p, and 


Pf were determined by extrapolation to 1,/L gR° =0. 
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Figure 1 - Shell Pressure Factor (p_) as a Function of Bulkhead Spacing (L,/R) 
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Figure 2 - Frame Pressure Factor (?) as a Function of Frame Stiffness (,/L ,R*) 
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E = + (100 h/R - | L,/k and 1,/L,R°. Here plots of constant values for p,, appear as fami- 


lies of curves intersecting along straight lines. These lines divide the graph into regions of 
n = 2, 3, 4, and 5. A close estimate of the collapse pressure can thus be obtained directly 
from this figure once the coordinates are known. For a more precise determination of p_ Fig- 


ures 1 and 2 should be used after the appropriate value of n has been found from Figure 3. 


ACCURACY OF RESULTS 


In the construction of Figures 1 through 3, some scatter in the calculated points was 
unavoidable, as might be expected. However, ‘in only a few isolated instances was the dif- 
ference between the graphical results and those from Kendrick Part If as high as 15 percent. 
In the great majority of cases, it was less than 10 percent. A representative set of results 


obtained from Figures 1 to 3 is presented for comparison with Kendrick’s theory in Table 1. 
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TABLE 1 


Comparison of Numerical Results with Kendrick’s (Part II) Theory 


: i Relative 
Relative Relative eranie 
Thickness | Length | ctitiness 
LAG). Teds 


4.93 x 10° 
5.99 x 10° 
6.62 x 10° 
7.15 x 10* 


4.93 x 10° 
5.99 x 10° 
6.62 x 10° 
7.15 x 107° 


4.93 x 10° 
5.99 x 10° 
6.62 x 10° 
7.15 x10° 


4.93 x 10° 
5.99 x 10° 
6.62 x 10° 
7.15 x10 


1.145 x 10° 
1.385 x 10° 
1.581 x 107° 
1.773 x 10-6 


1.145 x 10° 
1.385 x 10-6 
1.581 x 1076 
1.773 x 107° 


1.145 x 10-6 
1,385 x 10° 
1.581 x 1076 
1.773 x 10-6 


1.145 x 1076 
1.385 x 107 
1,581 x 197-6 
1.773 x 1076 


*For these three geometries, collapse pressures for n = 2 were nearly the same as 
forn = 3. 


Collapse Pressure p,, Mode of Collapse n 


From From From From 
Figures Kendrick’s Figure Kendrick’s 
1 and 2 Theory Theory 


BRO PO PO RM GPM PM WwW Ww S&S SP SBS PS POM MM MH NMP PY 


APPENDIX A 
NUMERICAL EXAMPLE 


A numerical example is provided to illustrate the use of the curves presented herein. 
A typical case has been chosen in which the dimensions could be those of a submarine pres- 


sure hull with external frames. The pertinent scantlings are: 


R = 96 in. Ar = 9.63 in.? 
h = 0.768 in. I, = 80.48 in.* 
Ly = 30 in. e = 4.49 in. 
L, = 384 in. 


The cylinder has 13 frames with one short bay at each end. First the value of the ef- 
fective length L. is found to be 13.48 in. Then 


I 
€_ is 6.85 x 107° 


3 
Lf 


L 
E ~ 1 (1004/R - vb is 4.26 
3 R 


From Figure 3, p,, is found to be approximately 1500 psi corresponding to the mode n = 3. 
Using the ratio L,/F = 4 in Figure 1, p, - R/100A is 150 psi. This is multiplied by 0.80, the 
value of 100A/R, to give p, = 120 psi. From Figure 2, the value Pr for n = 3 is 1370 psi. 


and 


Adding these two pressures, p_, is found to be 1490 psi as compared with 1500 psi found di- 
rectly from Figure 3. 

Table 2 compares the values of p, for this example as determined by Kendrick Part I, 
Kendrick Part III, Bryant, and this graphical method. While the pressure obtained from the 
graphs agrees closely with Kendrick Part III, it is interesting to note that both the Kendrick 


Part I and Bryant values are considerably higher and predict a different number of lobes. 


TABLE 2 


Example of Results Given by Various Methods 


Kendrick } Graphical | Kendrick 
Part | Method Part Ill 


EI 1855, 1490 1409 


2045 


Bryant 


APPENDIX B 
APPROXIMATE FORMULA FOR /, 


Considerable computation is required to determine the parameter I,, which is a measure 
of frame strength. To simplify this work, an approximate formula has been derived which is 
applicable to frames of various cross-sectional shapes. 

As previously defined, J, is the moment of inertia of the combination of a frame plus an 
effective length of shell L, about the centroid of the combined section. Figure 4 is a diagram 
of the most general frame shape in practical use, i.e., the builtup H-section with unequal 
flanges. Other shapes, such as T- or [-sections, are special cases of this shape. Referring 


to Figure 4, the quantity e of Equation [1] on page 8 can be written: 


1 be (« t ) as a ; cu2 h 
e=—— —— = = == —s) 
A; 5 + 5 (d+u yo + 5 + i [3] 


where a, 6, c, d, 8, t, and wu are defined in Figure 4. Using the relationship 


Ay = a8 + bt + cu [4] 


Equation [3] can be rearranged so that 
1 as uA; bt h 
face (d—?) Wiser. + —L + on (¢—w) + = [5] 


If the small quantity ¢ — u is neglected, Equation [5] becomes 


Be ane ey ancien Se h 
Ockre (d—t) 4 +uth. [6] 


2 
Aye? ~Z ao? [1 : 7 | 194 


and Ae? is written 


The quantity I; can also be simplified if a few small terms are neglected. The exact expres- 


sion for If is 


2 


1 ; 1 ‘ 
+ as | punerwam =e von | Sermyne| 


10 


2 
t h 
po (bt3 + a3s + cu?) + bt d-—--e+ — [8] 
12 2 


Using the expression for e of Equation [6] the above’ becomes 


1 bt (bt - cu) | ? 
= 3 3 
Vigaieees + a°s + cu? ) Nay peer | [9] 
és 2 19 2 
_ as Gat) (bt — cu) F cu (d— 02 1a bt -— cu 
4 Ay 4 Ar 


With the approximation wu = ¢, Equation [9] becomes 


1 bt bt - 2 
I, = — (bt? + as + cu3) + — (d—t)? 1- Z| 
12 4 


A; 


(bt — cu) cu bt —cu |? 
(d—¢) ——— i + — (dt)? [3 
a 4 ‘ 4; 


[10] 


or, 


1 a Ne 
po = (EP a eaw art) 4 SEZ" Fg ree lar) {11] 
12 4 


With a further simplification that the first term of Equation [11] is approximately 7 5 (a = He. 


the resulting expression for I is 


(d— t)? (bt + cw) 3 


po & a - a on? [12] 


In calculating /, from Equation [1], the term L,h3/12 is usually so small that it can be neg- 
lected. With the substitution of Equations [7] and [12] into Equation [1], the terms can be 


rearranged so that 


(d—t)? 1 bt — cu Atu |? 
soar 
i= A 1+ + 

jasc 


For an I-beam a simpler expression is obtained. Since 6 = c and u = ¢ Equation [3] 


becomes 
a+h 
e= 
9 [14] 
while the Equation [12] for I reduces to 
d — t)2 4bt 
ly = Css Ay 1 ees ({15] 
12 Ay 
The expression for /, is 
d+h)? d-t)? 4b¢ 
gel Grae ye) (1+ [16] 
4 Ay 3 
1+ Z 
L, 
For a T-section, the quantities c and wu are zero and Equation [5] can be written 
1 2 > 
eo Urey ete) 2 aa ees [17] 
A, ne 2 2 
which reduces to 
d bt t h 
eo (11> -—) + — [18] 
2 d 2 
Substituting this expression in Equation [8) 
1 da? bt 
I, = — (0¢3)+ a?s) + — ue (1-—-) [19] 
12 4 
or 
d? Qt bt t2 
fee 2) eee 
12 d A, da 
[20] 


12 


If the small quantity ¢?/d? is neglected, this can be reduced to 


2 
sept lia) (1-2 7) [21] 


i A, d 
and 
d2 b A-t \2 
ee (a: + (1-2) ( -F +s) [22] 
12 f d A, d Ay 


To demonstrate the accuracy of these 


frame formulas some numerical examples are 
provided. Table 3 lists the dimensions of 
several frames representing shapes and sizes 


which might be used in submarine construc- 


tion. Type 1 is an H-section with unequal 
flanges which was used as the example in 
Appendix A. Types 2 - 5 are I- and 
T-sections taken from the steel construction 


manual of the American Institute of Steel 


Construction.’ Table 4 lists shell dimen- 


sions for each of these frames and compares 


the values of I; and I calculated from the 


simplified formulas with the exact values from 
Reference 7. The error in each case is found Figure 4 - [Identification of Symbols 


to be less than 2 percent. Used in Calculating /, 


13 


TABLE 3 


Dimensions of Typical Frame Cross Sections 


ae. Frame Dimensions in inches 


as a 
Re ae 


Pree Perper f= Pero [= 
P sre [ aes | eco [soo [ — | osm | ose 


*Values of A, in inchese for Frames 2-5 are those given in the A.LS.C. Handbook. 2 In some 
cases, these include fillets not listed in this table. 


TABLE 4 


Comparison of Approximate and Exact Values for Ur and /, 


ky hy 
Shell Shell 
Radius | Thickness 
in in 


inches inches Approximate 


*Exact values of I, for Frames 2-5 are those given in Reference 7. 
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